Code
dtx <- growth %>%
mutate (strains = factor (strains
, levels = c ("Control" , "BR-11001" , "BR-11002" )
, ordered = TRUE )) %>%
merge (. , gsm) %>%
mutate (svi = (shoot_length + root_length)* grp)
# trait <- "shoot_length"
rs <- xfun:: cache_rds ({
traits <- names (dtx)[4 : length (dtx)]
results <- traits %>% set_names () %>% map ( \(trait) {
# Liberar memoria
gc ()
# Definir modelo
mdf <- paste (trait, "0 + strains*nacl + (1|repetition)" , sep = " ~ " )
# Mostrar qué rasgo se analiza
cat (" \n\n ## Análisis para: " , trait, " \n\n " )
cat (" \n\n ### Diagnostico de modelo \n\n " )
# Remover outliers (asumo que tienes esta función personalizada)
md <- dtx %>%
remove_outliers (as.formula (mdf),
plot_diag = TRUE ,
drop_na = TRUE )
md$ diagplot %>% print ()
if (nrow (md$ outliers) > 0 ) {
cat ("> Outliers detectados \n " )
md$ outliers %>% kable () %>% print ()
} else {
cat ("> No se detectaron outliers \n " )
}
mdc <- md$ model$ clean
cat (" \n\n ### ANOVA tipo III \n\n " )
Anova (mdc, type = 3 , test.statistic = "F" ) %>% print ()
cat (" \n\n ### Comparaciones múltiples (Emmeans) \n\n " )
# Comparaciones múltiples into the groups
mc_in <- emmeans (mdc, ~ strains| nacl, type = "response" ) %>%
cld (Letters = letters, reversed = TRUE ) %>%
mutate (across (".group" , trimws)) %>%
mutate (across (".group" , tolower)) %>%
rename (sig_in = ".group" )
mc_ou <- emmeans (mdc, ~ nacl| strains, type = "response" ) %>%
cld (Letters = letters, reversed = TRUE ) %>%
mutate (across (".group" , trimws)) %>%
mutate (across (".group" , toupper)) %>%
rename (sig_ou = ".group" )
mc <- merge (mc_in, mc_ou) %>%
unite (col = "sig" , c ("sig_in" , "sig_ou" ), sep = "" )
mc %>% kable () %>% print ()
cat (" \n\n ### Gráfico resumen \n\n " )
plot <- mc %>%
plot_smr (type = "line" ,
x = "nacl" ,
y = "emmean"
, group = "strains"
, ylab = trait
, sig = "sig"
, error = "SE" ) +
scale_color_manual (values = c ("#fbb4ae" , "#b3cde3" , "#ccebc5" ))
plot %>% print ()
list (mc = mc,
plot = plot
)
})
results
}, file = "result.rds" , rerun = T)
##
##
## ## Análisis para: shoot_length
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: shoot_length
## F Df Df.res Pr(>F)
## strains 51368.34 3 33.245 < 0.00000000000000022 ***
## nacl 6868.36 3 33.000 < 0.00000000000000022 ***
## strains:nacl 86.68 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--:|--------:|--------:|:---|
## |BR-11001 |0 | 5.809725| 0.0244621| 36| 5.760114| 5.859337|bA |
## |BR-11001 |150 | 5.797650| 0.0244621| 36| 5.748038| 5.847262|bA |
## |BR-11001 |300 | 4.558875| 0.0244621| 36| 4.509264| 4.608486|bB |
## |BR-11001 |450 | 3.062825| 0.0244621| 36| 3.013213| 3.112436|bC |
## |BR-11002 |0 | 6.020650| 0.0244621| 36| 5.971038| 6.070261|aA |
## |BR-11002 |150 | 6.033325| 0.0244621| 36| 5.983714| 6.082936|aA |
## |BR-11002 |300 | 5.262625| 0.0244621| 36| 5.213013| 5.312237|aB |
## |BR-11002 |450 | 3.959200| 0.0244621| 36| 3.909588| 4.008812|aC |
## |Control |0 | 4.770450| 0.0244621| 36| 4.720839| 4.820062|cA |
## |Control |150 | 4.727025| 0.0244621| 36| 4.677414| 4.776637|cA |
## |Control |300 | 4.136000| 0.0244621| 36| 4.086388| 4.185612|cB |
## |Control |450 | 2.110125| 0.0244621| 36| 2.060513| 2.159737|cC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: root_length
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: root_length
## F Df Df.res Pr(>F)
## strains 16443.532 3 31.751 < 0.00000000000000022 ***
## nacl 2541.001 3 33.000 < 0.00000000000000022 ***
## strains:nacl 86.177 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--------:|---------:|---------:|:---|
## |BR-11001 |0 | 3.316350| 0.0260871| 35.67545| 3.2634262| 3.3692738|bA |
## |BR-11001 |150 | 3.299150| 0.0260871| 35.67545| 3.2462262| 3.3520738|bA |
## |BR-11001 |300 | 2.493275| 0.0260871| 35.67545| 2.4403512| 2.5461988|bB |
## |BR-11001 |450 | 2.097775| 0.0260871| 35.67545| 2.0448512| 2.1506988|bC |
## |BR-11002 |0 | 3.934625| 0.0260871| 35.67545| 3.8817012| 3.9875488|aA |
## |BR-11002 |150 | 3.983000| 0.0260871| 35.67545| 3.9300762| 4.0359238|aA |
## |BR-11002 |300 | 3.457075| 0.0260871| 35.67545| 3.4041512| 3.5099988|aB |
## |BR-11002 |450 | 2.543650| 0.0260871| 35.67545| 2.4907262| 2.5965738|aC |
## |Control |0 | 3.051550| 0.0260871| 35.67545| 2.9986262| 3.1044738|cA |
## |Control |150 | 2.999075| 0.0260871| 35.67545| 2.9461512| 3.0519988|cA |
## |Control |300 | 2.153075| 0.0260871| 35.67545| 2.1001512| 2.2059988|cB |
## |Control |450 | 0.933900| 0.0260871| 35.67545| 0.8809762| 0.9868238|cC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: fresh_weight
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: fresh_weight
## F Df Df.res Pr(>F)
## strains 2441756 3 33.245 < 0.00000000000000022 ***
## nacl 337881 3 33.000 < 0.00000000000000022 ***
## strains:nacl 23722 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--:|---------:|---------:|:---|
## |BR-11001 |0 | 0.413850| 0.0002851| 36| 0.4132718| 0.4144282|bA |
## |BR-11001 |150 | 0.414225| 0.0002851| 36| 0.4136468| 0.4148032|bA |
## |BR-11001 |300 | 0.299875| 0.0002851| 36| 0.2992968| 0.3004532|bB |
## |BR-11001 |450 | 0.257300| 0.0002851| 36| 0.2567218| 0.2578782|bC |
## |BR-11002 |0 | 0.556675| 0.0002851| 36| 0.5560968| 0.5572532|aA |
## |BR-11002 |150 | 0.556750| 0.0002851| 36| 0.5561718| 0.5573282|aA |
## |BR-11002 |300 | 0.314775| 0.0002851| 36| 0.3141968| 0.3153532|aB |
## |BR-11002 |450 | 0.300075| 0.0002851| 36| 0.2994968| 0.3006532|aC |
## |Control |0 | 0.343000| 0.0002851| 36| 0.3424218| 0.3435782|cA |
## |Control |150 | 0.342850| 0.0002851| 36| 0.3422718| 0.3434282|cA |
## |Control |300 | 0.257125| 0.0002851| 36| 0.2565468| 0.2577032|cB |
## |Control |450 | 0.214125| 0.0002851| 36| 0.2135468| 0.2147032|cC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: dry_weight
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: dry_weight
## F Df Df.res Pr(>F)
## strains 14802.65 3 33.245 < 0.00000000000000022 ***
## nacl 1129.13 3 33.000 < 0.00000000000000022 ***
## strains:nacl 173.93 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--:|---------:|---------:|:---|
## |BR-11001 |0 | 0.032300| 0.0003048| 36| 0.0316818| 0.0329182|bA |
## |BR-11001 |150 | 0.032475| 0.0003048| 36| 0.0318568| 0.0330932|bA |
## |BR-11001 |300 | 0.027275| 0.0003048| 36| 0.0266568| 0.0278932|bB |
## |BR-11001 |450 | 0.025775| 0.0003048| 36| 0.0251568| 0.0263932|bC |
## |BR-11002 |0 | 0.049100| 0.0003048| 36| 0.0484818| 0.0497182|aA |
## |BR-11002 |150 | 0.049150| 0.0003048| 36| 0.0485318| 0.0497682|aA |
## |BR-11002 |300 | 0.030475| 0.0003048| 36| 0.0298568| 0.0310932|aB |
## |BR-11002 |450 | 0.031550| 0.0003048| 36| 0.0309318| 0.0321682|aB |
## |Control |0 | 0.026375| 0.0003048| 36| 0.0257568| 0.0269932|cA |
## |Control |150 | 0.026350| 0.0003048| 36| 0.0257318| 0.0269682|cA |
## |Control |300 | 0.020850| 0.0003048| 36| 0.0202318| 0.0214682|cB |
## |Control |450 | 0.018500| 0.0003048| 36| 0.0178818| 0.0191182|cC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: ms_porc
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: ms_porc
## F Df Df.res Pr(>F)
## strains 7712.487 3 32.162 < 0.00000000000000022 ***
## nacl 226.923 3 33.000 < 0.00000000000000022 ***
## strains:nacl 12.429 6 33.000 0.0000002787 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|---------:|---------:|--------:|---------:|---------:|:---|
## |BR-11001 |0 | 7.804621| 0.0901902| 35.81973| 7.621675| 7.987568|bC |
## |BR-11001 |150 | 7.839800| 0.0901902| 35.81973| 7.656854| 8.022746|bC |
## |BR-11001 |300 | 9.094932| 0.0901902| 35.81973| 8.911986| 9.277878|bB |
## |BR-11001 |450 | 10.017337| 0.0901902| 35.81973| 9.834391| 10.200283|bA |
## |BR-11002 |0 | 8.820208| 0.0901902| 35.81973| 8.637262| 9.003155|aC |
## |BR-11002 |150 | 8.827945| 0.0901902| 35.81973| 8.644999| 9.010891|aC |
## |BR-11002 |300 | 9.681266| 0.0901902| 35.81973| 9.498320| 9.864212|aB |
## |BR-11002 |450 | 10.513656| 0.0901902| 35.81973| 10.330710| 10.696602|aA |
## |Control |0 | 7.689417| 0.0901902| 35.81973| 7.506471| 7.872364|bC |
## |Control |150 | 7.685450| 0.0901902| 35.81973| 7.502504| 7.868397|bC |
## |Control |300 | 8.108541| 0.0901902| 35.81973| 7.925594| 8.291487|cB |
## |Control |450 | 8.639402| 0.0901902| 35.81973| 8.456456| 8.822348|cA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: humidity_cont
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: humidity_cont
## F Df Df.res Pr(>F)
## strains 11634347 3 33.245 < 0.00000000000000022 ***
## nacl 1686160 3 33.000 < 0.00000000000000022 ***
## strains:nacl 111712 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--:|---------:|---------:|:---|
## |BR-11001 |0 | 0.381550| 0.0001198| 36| 0.3813071| 0.3817929|bA |
## |BR-11001 |150 | 0.381750| 0.0001198| 36| 0.3815071| 0.3819929|bA |
## |BR-11001 |300 | 0.272600| 0.0001198| 36| 0.2723571| 0.2728429|bB |
## |BR-11001 |450 | 0.231525| 0.0001198| 36| 0.2312821| 0.2317679|bC |
## |BR-11002 |0 | 0.507575| 0.0001198| 36| 0.5073321| 0.5078179|aA |
## |BR-11002 |150 | 0.507600| 0.0001198| 36| 0.5073571| 0.5078429|aA |
## |BR-11002 |300 | 0.284300| 0.0001198| 36| 0.2840571| 0.2845429|aB |
## |BR-11002 |450 | 0.268525| 0.0001198| 36| 0.2682821| 0.2687679|aC |
## |Control |0 | 0.316625| 0.0001198| 36| 0.3163821| 0.3168679|cA |
## |Control |150 | 0.316500| 0.0001198| 36| 0.3162571| 0.3167429|cA |
## |Control |300 | 0.236275| 0.0001198| 36| 0.2360321| 0.2365179|cB |
## |Control |450 | 0.195625| 0.0001198| 36| 0.1953821| 0.1958679|cC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: humidity_porc
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: humidity_porc
## F Df Df.res Pr(>F)
## strains 987417.131 3 32.162 < 0.00000000000000022 ***
## nacl 226.923 3 33.000 < 0.00000000000000022 ***
## strains:nacl 12.429 6 33.000 0.0000002787 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--------:|--------:|--------:|:---|
## |BR-11001 |0 | 92.19538| 0.0901902| 35.81973| 92.01243| 92.37832|aA |
## |BR-11001 |150 | 92.16020| 0.0901902| 35.81973| 91.97725| 92.34315|aA |
## |BR-11001 |300 | 90.90507| 0.0901902| 35.81973| 90.72212| 91.08801|bB |
## |BR-11001 |450 | 89.98266| 0.0901902| 35.81973| 89.79972| 90.16561|bC |
## |BR-11002 |0 | 91.17979| 0.0901902| 35.81973| 90.99685| 91.36274|bA |
## |BR-11002 |150 | 91.17205| 0.0901902| 35.81973| 90.98911| 91.35500|bA |
## |BR-11002 |300 | 90.31873| 0.0901902| 35.81973| 90.13579| 90.50168|cB |
## |BR-11002 |450 | 89.48634| 0.0901902| 35.81973| 89.30340| 89.66929|cC |
## |Control |0 | 92.31058| 0.0901902| 35.81973| 92.12764| 92.49353|aA |
## |Control |150 | 92.31455| 0.0901902| 35.81973| 92.13160| 92.49750|aA |
## |Control |300 | 91.89146| 0.0901902| 35.81973| 91.70851| 92.07441|aB |
## |Control |450 | 91.36060| 0.0901902| 35.81973| 91.17765| 91.54354|aC |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: sod
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: sod
## F Df Df.res Pr(>F)
## strains 869.82 3 11.003 0.0000000000002362 ***
## nacl 57039.15 3 33.000 < 0.00000000000000022 ***
## strains:nacl 310.48 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--------:|--------:|--------:|:---|
## |BR-11001 |0 | 26.78000| 0.5248219| 3.291193| 25.19035| 28.36965|aD |
## |BR-11001 |150 | 39.41000| 0.5248219| 3.291193| 37.82035| 40.99965|bC |
## |BR-11001 |300 | 53.43625| 0.5248219| 3.291193| 51.84660| 55.02590|bB |
## |BR-11001 |450 | 63.35875| 0.5248219| 3.291193| 61.76910| 64.94840|bA |
## |BR-11002 |0 | 26.67125| 0.5248219| 3.291193| 25.08160| 28.26090|aD |
## |BR-11002 |150 | 42.10250| 0.5248219| 3.291193| 40.51285| 43.69215|aC |
## |BR-11002 |300 | 57.90625| 0.5248219| 3.291193| 56.31660| 59.49590|aB |
## |BR-11002 |450 | 68.11625| 0.5248219| 3.291193| 66.52660| 69.70590|aA |
## |Control |0 | 23.96375| 0.5248219| 3.291193| 22.37410| 25.55340|bD |
## |Control |150 | 36.23125| 0.5248219| 3.291193| 34.64160| 37.82090|cC |
## |Control |300 | 48.24000| 0.5248219| 3.291193| 46.65035| 49.82965|cB |
## |Control |450 | 56.25375| 0.5248219| 3.291193| 54.66410| 57.84340|cA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: cat
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: cat
## F Df Df.res Pr(>F)
## strains 392.73 3 11.136 0.00000000001432 ***
## nacl 26515.18 3 33.000 < 0.00000000000000022 ***
## strains:nacl 231.98 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--------:|--------:|--------:|:---|
## |BR-11001 |0 | 14.93625| 0.4490014| 3.456453| 13.60837| 16.26413|bD |
## |BR-11001 |150 | 25.46625| 0.4490014| 3.456453| 24.13837| 26.79413|bC |
## |BR-11001 |300 | 35.49500| 0.4490014| 3.456453| 34.16712| 36.82288|bB |
## |BR-11001 |450 | 41.70000| 0.4490014| 3.456453| 40.37212| 43.02788|bA |
## |BR-11002 |0 | 16.01250| 0.4490014| 3.456453| 14.68462| 17.34038|aD |
## |BR-11002 |150 | 26.80125| 0.4490014| 3.456453| 25.47337| 28.12913|aC |
## |BR-11002 |300 | 39.09125| 0.4490014| 3.456453| 37.76337| 40.41913|aB |
## |BR-11002 |450 | 46.05875| 0.4490014| 3.456453| 44.73087| 47.38663|aA |
## |Control |0 | 14.15125| 0.4490014| 3.456453| 12.82337| 15.47913|cD |
## |Control |150 | 21.94750| 0.4490014| 3.456453| 20.61962| 23.27538|cC |
## |Control |300 | 31.12375| 0.4490014| 3.456453| 29.79587| 32.45163|cB |
## |Control |450 | 35.72500| 0.4490014| 3.456453| 34.39712| 37.05288|cA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: apx
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: apx
## F Df Df.res Pr(>F)
## strains 142.535 3 11.265 0.000000003101 ***
## nacl 8002.578 3 33.000 < 0.00000000000000022 ***
## strains:nacl 70.101 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--------:|---------:|--------:|:---|
## |BR-11001 |0 | 10.84000| 0.5221022| 3.619693| 9.328325| 12.35167|aD |
## |BR-11001 |150 | 17.74625| 0.5221022| 3.619693| 16.234575| 19.25792|bC |
## |BR-11001 |300 | 25.33125| 0.5221022| 3.619693| 23.819575| 26.84292|bB |
## |BR-11001 |450 | 29.64125| 0.5221022| 3.619693| 28.129575| 31.15292|bA |
## |BR-11002 |0 | 10.75625| 0.5221022| 3.619693| 9.244575| 12.26792|aD |
## |BR-11002 |150 | 19.40750| 0.5221022| 3.619693| 17.895825| 20.91917|aC |
## |BR-11002 |300 | 28.38000| 0.5221022| 3.619693| 26.868325| 29.89167|aB |
## |BR-11002 |450 | 33.04625| 0.5221022| 3.619693| 31.534575| 34.55792|aA |
## |Control |0 | 9.36000| 0.5221022| 3.619693| 7.848325| 10.87167|bD |
## |Control |150 | 15.19750| 0.5221022| 3.619693| 13.685825| 16.70917|cC |
## |Control |300 | 21.74500| 0.5221022| 3.619693| 20.233325| 23.25667|cB |
## |Control |450 | 25.58750| 0.5221022| 3.619693| 24.075825| 27.09917|cA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: gpx
##
##
##
## ### Diagnostico de modelo
## > Outliers detectados
##
##
## | |index |strains |nacl |repetition | gpx| resi| res_MAD| rawp.BHStud| adjp| bholm|out_flag |
## |:--|:-----|:--------|:----|:----------|----:|----------:|---------:|-----------:|------------:|--------:|:--------|
## |20 |20 |BR-11002 |0 |R4 | 7.74| -0.7986443| -3.515746| 0.0004385| 0.0004385205| 0.021049|OUTLIER |
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: gpx
## F Df Df.res Pr(>F)
## strains 65.977 3 10.985 0.0000002597 ***
## nacl 5923.719 3 32.001 < 0.00000000000000022 ***
## strains:nacl 69.612 6 32.001 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|---------:|---------:|--------:|---------:|---------:|:---|
## |BR-11001 |0 | 7.933750| 0.5459773| 3.319524| 6.287096| 9.580404|aD |
## |BR-11001 |150 | 11.718750| 0.5459773| 3.319524| 10.072096| 13.365404|bC |
## |BR-11001 |300 | 17.838750| 0.5459773| 3.319524| 16.192096| 19.485404|bB |
## |BR-11001 |450 | 20.861250| 0.5459773| 3.319524| 19.214596| 22.507904|bA |
## |BR-11002 |0 | 7.615149| 0.5513238| 3.450493| 5.983326| 9.246973|abD |
## |BR-11002 |150 | 13.461250| 0.5459773| 3.319524| 11.814596| 15.107904|aC |
## |BR-11002 |300 | 18.953750| 0.5459773| 3.319524| 17.307096| 20.600404|aB |
## |BR-11002 |450 | 23.381250| 0.5459773| 3.319524| 21.734596| 25.027904|aA |
## |Control |0 | 7.148750| 0.5459773| 3.319524| 5.502096| 8.795404|bD |
## |Control |150 | 10.266250| 0.5459773| 3.319524| 8.619596| 11.912904|cC |
## |Control |300 | 14.838750| 0.5459773| 3.319524| 13.192096| 16.485404|cB |
## |Control |450 | 17.778750| 0.5459773| 3.319524| 16.132096| 19.425404|cA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: grp
##
##
##
## ### Diagnostico de modelo
## > Outliers detectados
##
##
## | |index |strains |nacl |repetition | grp| resi| res_MAD| rawp.BHStud| adjp| bholm|out_flag |
## |:--|:-----|:--------|:----|:----------|---:|----:|----------------:|-----------:|----:|-----:|:--------|
## |9 |9 |BR-11001 |300 |R1 | 98| -0.5| -23731533851904| 0| 0| 0|OUTLIER |
## |10 |10 |BR-11001 |300 |R2 | 100| 1.5| 71194601555710| 0| 0| 0|OUTLIER |
## |11 |11 |BR-11001 |300 |R3 | 98| -0.5| -23731533851904| 0| 0| 0|OUTLIER |
## |12 |12 |BR-11001 |300 |R4 | 98| -0.5| -23731533851904| 0| 0| 0|OUTLIER |
## |13 |13 |BR-11001 |450 |R1 | 76| -2.5| -118657669259517| 0| 0| 0|OUTLIER |
## |14 |14 |BR-11001 |450 |R2 | 78| -0.5| -23731533851903| 0| 0| 0|OUTLIER |
## |15 |15 |BR-11001 |450 |R3 | 80| 1.5| 71194601555710| 0| 0| 0|OUTLIER |
## |16 |16 |BR-11001 |450 |R4 | 80| 1.5| 71194601555710| 0| 0| 0|OUTLIER |
## |25 |25 |BR-11002 |300 |R1 | 98| -1.0| -47463067703807| 0| 0| 0|OUTLIER |
## |26 |26 |BR-11002 |300 |R2 | 98| -1.0| -47463067703807| 0| 0| 0|OUTLIER |
## |27 |27 |BR-11002 |300 |R3 | 100| 1.0| 47463067703807| 0| 0| 0|OUTLIER |
## |28 |28 |BR-11002 |300 |R4 | 100| 1.0| 47463067703807| 0| 0| 0|OUTLIER |
## |29 |29 |BR-11002 |450 |R1 | 82| -2.0| -94926135407613| 0| 0| 0|OUTLIER |
## |31 |31 |BR-11002 |450 |R3 | 86| 2.0| 94926135407614| 0| 0| 0|OUTLIER |
## |41 |41 |Control |300 |R1 | 96| 0.5| 23731533851903| 0| 0| 0|OUTLIER |
## |42 |42 |Control |300 |R2 | 94| -1.5| -71194601555711| 0| 0| 0|OUTLIER |
## |43 |43 |Control |300 |R3 | 96| 0.5| 23731533851903| 0| 0| 0|OUTLIER |
## |44 |44 |Control |300 |R4 | 96| 0.5| 23731533851903| 0| 0| 0|OUTLIER |
## |45 |45 |Control |450 |R1 | 56| 3.5| 166120736963325| 0| 0| 0|OUTLIER |
## |46 |46 |Control |450 |R2 | 52| -0.5| -23731533851903| 0| 0| 0|OUTLIER |
## |47 |47 |Control |450 |R3 | 52| -0.5| -23731533851903| 0| 0| 0|OUTLIER |
## |48 |48 |Control |450 |R4 | 50| -2.5| -118657669259517| 0| 0| 0|OUTLIER |
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: grp
## F Df Df.res
## strains 89732522025707554560428464086064.0000 3 13.111
## nacl 674537293467640588642442006462.0000 2 16.361
## strains:nacl 2.5243 2 16.003
## Pr(>F)
## strains <0.0000000000000002 ***
## nacl <0.0000000000000002 ***
## strains:nacl 0.1115
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|------:|--:|--------:|--------:|--------:|:---|
## |BR-11001 |0 | 100| 0| 15.95945| 100| 100|aA |
## |BR-11001 |150 | 100| 0| 15.95945| 100| 100|aA |
## |BR-11001 |450 | | | | | | |
## |BR-11002 |0 | 100| 0| 15.95945| 100| 100|aA |
## |BR-11002 |150 | 100| 0| 15.95945| 100| 100|aA |
## |BR-11002 |450 | 84| 0| 18.78373| 84| 84|aB |
## |Control |0 | 100| 0| 15.95945| 100| 100|aA |
## |Control |150 | 100| 0| 15.95945| 100| 100|aA |
## |Control |450 | | | | | | |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: mgt
##
##
##
## ### Diagnostico de modelo
## > Outliers detectados
##
##
## | |index |strains |nacl |repetition | mgt| resi| res_MAD| rawp.BHStud| adjp| bholm|out_flag |
## |:--|:-----|:--------|:----|:----------|--------:|----------:|---------:|-----------:|--------------:|---------:|:--------|
## |30 |30 |BR-11002 |450 |R2 | 3.166667| -0.4404369| -4.816386| 0.0000015| 0.000001461813| 0.0000702|OUTLIER |
## |47 |47 |Control |450 |R3 | 6.115385| 0.3376700| 3.692582| 0.0002220| 0.000221988925| 0.0104335|OUTLIER |
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: mgt
## F Df Df.res Pr(>F)
## strains 592.599 3 31.233 < 0.00000000000000022 ***
## nacl 2874.590 3 31.176 < 0.00000000000000022 ***
## strains:nacl 50.111 6 31.279 0.00000000000001103 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|--------:|---------:|--:|--------:|--------:|:---|
## |BR-11001 |0 | 1.150000| 0.0495327| 34| 1.049338| 1.250663|bD |
## |BR-11001 |150 | 1.700000| 0.0495327| 34| 1.599338| 1.800663|aC |
## |BR-11001 |300 | 2.125816| 0.0495327| 34| 2.025154| 2.226479|bB |
## |BR-11001 |450 | 5.216262| 0.0495327| 34| 5.115599| 5.316924|bA |
## |BR-11002 |0 | 1.110000| 0.0495327| 34| 1.009337| 1.210662|bD |
## |BR-11002 |150 | 1.405000| 0.0495327| 34| 1.304337| 1.505662|bC |
## |BR-11002 |300 | 1.741429| 0.0495327| 34| 1.640766| 1.842091|cB |
## |BR-11002 |450 | 3.801322| 0.0580728| 34| 3.683304| 3.919340|cA |
## |Control |0 | 1.355000| 0.0495327| 34| 1.254338| 1.455663|aD |
## |Control |150 | 1.750000| 0.0495327| 34| 1.649337| 1.850663|aC |
## |Control |300 | 2.459774| 0.0495327| 34| 2.359111| 2.560436|aB |
## |Control |450 | 5.617582| 0.0580728| 34| 5.499564| 5.735601|aA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: unc
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: unc
## F Df Df.res Pr(>F)
## strains 65.2786 3 23.393 0.00000000001624 ***
## nacl 806.2880 3 33.000 < 0.00000000000000022 ***
## strains:nacl 2.9991 6 33.000 0.01887 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|---------:|---------:|--------:|---------:|---------:|:---|
## |BR-11001 |0 | 0.5316838| 0.0649248| 26.43367| 0.3983355| 0.6650321|bD |
## |BR-11001 |150 | 1.2286803| 0.0649248| 26.43367| 1.0953321| 1.3620286|aC |
## |BR-11001 |300 | 1.6943196| 0.0649248| 26.43367| 1.5609713| 1.8276679|aB |
## |BR-11001 |450 | 2.9065733| 0.0649248| 26.43367| 2.7732250| 3.0399216|aA |
## |BR-11002 |0 | 0.4961936| 0.0649248| 26.43367| 0.3628453| 0.6295419|bD |
## |BR-11002 |150 | 0.9185542| 0.0649248| 26.43367| 0.7852059| 1.0519025|bC |
## |BR-11002 |300 | 1.5535492| 0.0649248| 26.43367| 1.4202010| 1.6868975|aB |
## |BR-11002 |450 | 2.8281305| 0.0649248| 26.43367| 2.6947822| 2.9614788|aA |
## |Control |0 | 0.7738005| 0.0649248| 26.43367| 0.6404522| 0.9071488|aD |
## |Control |150 | 1.3523896| 0.0649248| 26.43367| 1.2190414| 1.4857379|aC |
## |Control |300 | 1.6980144| 0.0649248| 26.43367| 1.5646661| 1.8313627|aB |
## |Control |450 | 2.8561455| 0.0649248| 26.43367| 2.7227973| 2.9894938|aA |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: syn
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: syn
## F Df Df.res Pr(>F)
## strains 1743.0005 3 23.735 < 0.00000000000000022 ***
## nacl 1196.5978 3 33.000 < 0.00000000000000022 ***
## strains:nacl 6.1427 6 33.000 0.0002109 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|---------:|---------:|--------:|---------:|---------:|:---|
## |BR-11001 |0 | 0.8293878| 0.0161694| 27.01915| 0.7962120| 0.8625635|aA |
## |BR-11001 |150 | 0.6267347| 0.0161694| 27.01915| 0.5935589| 0.6599105|bB |
## |BR-11001 |300 | 0.3208248| 0.0161694| 27.01915| 0.2876491| 0.3540006|bC |
## |BR-11001 |450 | 0.1315092| 0.0161694| 27.01915| 0.0983334| 0.1646850|aD |
## |BR-11002 |0 | 0.8012245| 0.0161694| 27.01915| 0.7680487| 0.8344003|aA |
## |BR-11002 |150 | 0.7157143| 0.0161694| 27.01915| 0.6825385| 0.7488901|aB |
## |BR-11002 |300 | 0.4118282| 0.0161694| 27.01915| 0.3786525| 0.4450040|aC |
## |BR-11002 |450 | 0.1491333| 0.0161694| 27.01915| 0.1159575| 0.1823091|aD |
## |Control |0 | 0.7502041| 0.0161694| 27.01915| 0.7170283| 0.7833799|bA |
## |Control |150 | 0.5695918| 0.0161694| 27.01915| 0.5364161| 0.6027676|cB |
## |Control |300 | 0.3441644| 0.0161694| 27.01915| 0.3109886| 0.3773401|bC |
## |Control |450 | 0.1170726| 0.0161694| 27.01915| 0.0838969| 0.1502484|aD |
##
##
## ### Gráfico resumen
##
##
## ## Análisis para: svi
##
##
##
## ### Diagnostico de modelo
## > No se detectaron outliers
##
##
## ### ANOVA tipo III
##
## Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
##
## Response: svi
## F Df Df.res Pr(>F)
## strains 31457.43 3 31.652 < 0.00000000000000022 ***
## nacl 7867.98 3 33.000 < 0.00000000000000022 ***
## strains:nacl 85.97 6 33.000 < 0.00000000000000022 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## ### Comparaciones múltiples (Emmeans)
##
##
##
## |strains |nacl | emmean| SE| df| lower.CL| upper.CL|sig |
## |:--------|:----|---------:|--------:|--------:|--------:|---------:|:---|
## |BR-11001 |0 | 912.6075| 4.911396| 35.63542| 902.6432| 922.5718|bA |
## |BR-11001 |150 | 909.6800| 4.911396| 35.63542| 899.7157| 919.6443|bA |
## |BR-11001 |300 | 694.6879| 4.911396| 35.63542| 684.7235| 704.6522|bB |
## |BR-11001 |450 | 405.1010| 4.911396| 35.63542| 395.1367| 415.0654|bC |
## |BR-11002 |0 | 995.5275| 4.911396| 35.63542| 985.5632| 1005.4918|aA |
## |BR-11002 |150 | 1001.6325| 4.911396| 35.63542| 991.6682| 1011.5968|aA |
## |BR-11002 |300 | 863.2957| 4.911396| 35.63542| 853.3314| 873.2600|aB |
## |BR-11002 |450 | 546.1985| 4.911396| 35.63542| 536.2342| 556.1629|aC |
## |Control |0 | 782.2000| 4.911396| 35.63542| 772.2357| 792.1643|cA |
## |Control |150 | 772.6100| 4.911396| 35.63542| 762.6457| 782.5743|cA |
## |Control |300 | 600.6111| 4.911396| 35.63542| 590.6468| 610.5754|cB |
## |Control |450 | 159.7628| 4.911396| 35.63542| 149.7984| 169.7271|cC |
##
##
## ### Gráfico resumen